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Abstract 

We study a percolation problem on a substrate formed by two-dimensional XY spin config- 
urations, using Monte Carlo methods. For a given spin configuration we construct percolation 
clusters by randomly choosing a direction x in the spin vector space, and then placing a perco- 
lation bond between nearest-neighbor sites i and j with probability pij = max(0, 1 — e * j ), 
where K > governs the percolation process. A line of percolation thresholds K C (J) is found 
in the low-temperature range J > J c , where J > is the XY coupling strength. Analysis of the 
correlation function g p (r), defined as the probability that two sites separated by a distance r belong 
to the same percolation cluster, yields algebraic decay for K > K C (J), and the associated critical 
exponent depends on J and K. Along the threshold line K C (J), the scaling dimension for g p is, 
within numerical uncertainties, equal to 1/8. On this basis, we conjecture that the percolation 
transition along the K C (J) line is of the Berezinskii-Kosterlitz-Thouless type. 

PACS numbers: 05.50.+q(lattice theory and statistics), 64.60.ah(percolation), 64.60.F-(equilibrium proper- 
ties near critical points, critical exponents), 75.10.Hk(classical spin models) 
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I. INTRODUCTION 



The XY model is formulated in terms of two-dimensional spins s normalized as \s\ — 1, 
residing on the sites of a lattice. The reduced Hamiltonian of the XY model (already divided 
by k-^T with k-Q the Boltzmann constant and T the temperature) reads 



H 



Si ■ s 



(1) 



where the sum is over all nearest-neighbor pairs, and J > is the ferromagnetic coupling 
strength. The spins are labeled by their site numbers. 

It is known from the Mermin-Wagner-Hohenberg-Coleman theorem [l| that there cannot 
exist spontaneous long-range order as long as J is finite in Eq. (CQ), because thermal fluctu- 
ations are strong enough to destroy the order. Nevertheless, the system undergoes a phase 
transition js — 4 1 as the coupling strength J increases. This type of transition is of infinite 
order and is known as the Berezinskii-Kosterlitz-Thouless (BKT) transition. For J < J c , the 
spin-spin correlation function decays exponentially, and the spins form a plasma of vortices; 
but for J > J c , the spin-spin correlation function decays algebraically with an exponent 
depending on J, and the spin configurations contain bound vortex-antivortex pairs. Tran- 
sitions of the BKT type occur in various kinds of systems. The XY-type of transition is 
related by duality to roughening transitions in solid-on-solid and related models J5|. Apart 

from the XY model, BKT transitions are found, among others, in vertex models (6|, models 

n fl ft 

of crystal surfaces [2J , the antiferromagnetic triangular Ising model [8] , string theory [9] , net- 



work systems jlO|, superfluid systems [ill ], and superconducting systems [l2|. These models 
may involve long-range or short-range interactions. 

It is also known that certain observables of statistical models are equivalent or closely 
related to properly defined geometric quantities. For instance, the Potts model can be 
exactly mapped onto the random-cluster model [lsj], and the susceptibility x °f the former 
is related to the cluster-size distribution of the latter; a similar situation applies to the 
Nienhuis O(n) loop model and the equivalent spin model [15]. The Mott-to-superfluid 



transition in the Bose-Hubbard model can be characterized by the winding number of the 



world lines of the particles [l6|] . The geometric percolation [l^] process has been employed to 



study percolation on critical substrates, such as the Ising model 
the O(n) model 22[ and even quantum Hall systems 23] . 
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-|20|, the Potts model [2l|, 
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In this work, we study the percolation problem on the substrate of the XY model ([[]). 
There is still some freedom in the choice of the percolation criterion. For instance, one may 
place percolation bonds between all neighboring XY spins if their orientations differ less 
than a given angle called the "conducting angle" . This problem was recently investigated by 
Wang et al. [24| . Here we use a different criterion. For a given spin configuration, we choose 
a randomly oriented Cartesian reference frame (x, y) in the two-dimensional spin space, and 
place bonds between nearest-neighbor pairs, say sites i and j, with a probability 

p 4 , =max(0,l-e- 2i ^), (2) 

where K > parametrizes the percolation problem. Note that for K = J these percola- 
tion clusters reduce to those formed by the cluster simulation process of the XY model as 
described in Sec. HT1 

The rest of this work is organized as follows. Section II III presents our numerical results 
for the critical points of the XY model on the square as well as on the triangular lattice. 
Section IIVI describes the analysis of the percolation problem for both lattices, with an em- 
phasis on the determination of the universal character of this type of percolation transition. 
We conclude with a discussion in Sec. [V] 



II. ALGORITHM AND SAMPLED QUANTITIES 
A. Spin updating algorithm 



We emp 



method 
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oy efficient Monte Carlo simulations of the XY model ([T]) by means of a cluster 



261 ] . We use a full cluster decomposition |25[ as follows: 



1. Choose a randomly oriented Cartesian frame of reference (x, y) in the spin space, and 
project the spin along the x and y axes as s = s x x + s y y. Accordingly, the scalar 
product in Eq. ([1]) is written s*j • Sj = s^s* + sfs V j, so that the Hamiltonian separates 
into two parts as H = W x + T-L y , with T-L x = — J sf and similar for % y . 

m 

2. Between each pair of nearest-neighboring sites, e.g., site i and j, place a bond with 
probability ptj = max(0, 1 — e~ 2JSi s i). 

3. Construct clusters on the basis of the occupied bonds. 



4. Independently for each cluster, flip the x components of all the spins in the cluster 
with probability 1/2. 

B. Sampled quantities 

We sampled several quantities, including the second and the fourth moments of the mag- 
netization density, M 2 = | ^uSk^jV 2 an d M A = | Ylk^l^/V 4 - These quantities determine 
the dimensionless Binder ratio [27| as 

Q m = (M 2 ) 2 /(M 4 ) , (3) 

where V = L 2 is the volume. The susceptibility is x = V(M 2 ). 

Denoting the size of the zth cluster by Cj, we also sampled the second and the fourth 
moments of the cluster size distribution as 

S2 = ^H c i and fif 4 = ^E c ?- ( 4 ) 

i i 

Accordingly, we define another dimensionless ratio as 

Q l = (S 2 ) 2 /(3(S 2 2 )-2(S 4 )). (5) 

Note that for the Ising model Qi in Eq. dSJ) is equal to Q m in Eq. (J3]). 

Also the spin-spin correlation g s (r) over distances r = L/2 and L/A was sampled. A third 
dimensionless ratio Q s is defined as 

Q s = (g s {L/2))/(g s {L/4)}. (6) 

In the high-temperature range J < J c , the spin-spin correlation decays exponentially, and Q s 
goes to as L — > oo. At criticality, however, g s {r) tends to algebraic decay as r~ 2xh , and Q s 
converges to a nontrivial universal value. In an ordered state with a non-zero magnetization 
density, Q s would converge to 1 instead. 

Finally, we define the correlation function g p (r) as the probability that two sites at a 
distance r belong to the same cluster. We sampled g p over distances r = L/2 and r = L/A. 
The associated dimensionless ratio is defined as 

Q p =(g p (L/2))/(g p (L/A)). (7) 
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III. CRITICAL POINTS 



A. Square lattice 

We simulated the XY model on L x L square lattices with periodic boundary conditions, 
with system sizes in the range 4 < L < 1024. As usual in Monte Carlo studies, the location 
of a critical point can well be determined using a dimensionless ratio. This is shown in 
Fig. [T]for the Binder ratio Q m . For J < J c , Q m approaches the infinite-temperature value 
1/2 as L — )■ oo, as expected for a normal distribution of the x and y components of the 
magnetization. For J > J c , Q m rapidly converges to a temperature-dependent value, as 
expected in the low-temperature XY phase. 




FIG. 1: (color online). Binder ratio Q m vs. coupling strength J for the square lattice. The lines 
connecting data points are added for clarity. 



Making use of the known magnetic scaling dimension Xh = 1/8 4j at the BKT transition, 
and the logarithmic correction factor with exponent 1/8 4j, |2j|, we expect that the scaled 
quantity xL 2Xh_2 (lnL) -1 / 8 tends to a constant at the transition point. The intersections in 
Fig. [21 which shows this scaled quantity as a function of J for several system sizes, confirm 
this expectation. 

Using the least-squares criterion, we fitted the quantity xL 2x ' i_2 (lnL) _1 / 8 data by the 
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FIG. 2: (color online). Scaled susceptibility xL 2zh_2 (lnL) -1 / 8 vs. coupling strength J for the 
square lattice, with Xh = 1/8. The lower figure is an enlarged version, and includes data for 
L = 256 and 512. The lines connecting the data points are added for clarity. 



formula 

3 2 

X L 2 ^~ 2 (lnL)-^ 8 = a + J2<Jc-JY(lnLY + J2rj(Jc-J) j 

i=i j=i 

+ ^/(lnLj+M^ + M,- 2 , (8) 

where the multiplicative and additive logarithmic corrections have been taken into account. 
We find that the data for 16 < L < 1024 and 1.100 < J < 1.125 are well described by 
Eq. ©. The fit yields J c = 1.124 (3). 

For the Ising model, one can prove that x = V(M 2 ) = V(S2), which exactly relates the 
thermodynamic quantity % to the geometric quantity S2. We thus expect that, in the case of 
the XY model, the singularity of £2 coincides with that of x- The data for S , 2 L 2:r ' i (lnL) _1 / 8 
(shown in Fig. [3]) were fitted by Eq. flH]). This fit yields J c = 1.120 (9), consistent with the 
result from x- 
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FIG. 3: (color online). Scaled second moment S2L 2xh (lnL)~ 1 / 8 of the cluster size distribution 
versus coupling strength J. These results apply to the cluster decomposition of the spin model on 
the square lattice. The lines connecting the data points are added for clarity. 



There exist already many estimates for the critical point of the XY model on the square 
lattice, the latest of which are J c = 1.1199(1) by Martin and Klaus |29j, J c = 1.1198 (14) by 
Butera and Pernici [30} , and J c = 1.1200 (1) by Arisue [3JJ. Our result for J c is consistent 
with these existing values. 



B. Triangular lattice 

We also simulated the XY model on the triangular lattice with periodic boundary con- 
ditions, for linear system sizes L in the range 4 < L < 512. The BKT phase transition 
is clearly exposed by Fig. HI which plots the ratio Q s = g s (r = L/2)/g s (r = L/4) versus 
the coupling strength J. In the high-temperature range J < J c , Q s rapidly approaches 
zero, which reflects the absence of long-range correlations; in the low-temperature range 
J > J c , it converges to a J-dependent value smaller than 1, in agreement with the presence 
of algebraically decaying correlations and the absence of a spontaneous magnetization. 

The g s (L/2)L 2xh (hiL)~' 1 / 8 data near criticality are shown in Fig. |5j They were fitted by 
Eq. dSJ), which yielded J c = 0.6833 (6). Analogous analyses were performed for the scaled 
susceptibility xL 2:E ' i_2 (lnL)~ 1 / 8 , leading to J c = 0.6831 (6). Our results for the critical 
coupling are consistent with the latest result J c = 0.6824 (8) by Butera and Pernici jsoj]. 
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J 

FIG. 4: (color online). Ratio Q s = g s (r = L/2)/g s (r = L/4) vs. coupling strength J for the 
triangular lattice, with g s the spin-spin correlation. The lines connecting data points are added for 
clarity. 

0.9 

0.85 

so 0.8 

^ 0.75 
-S 

§~ °- 7 
0.65 

0.6 

0.55 

0.66 0.68 0.7 0.72 0.74 

J 

FIG. 5: (color online). Scaled correlation g s (r = L/2)L 2xh (In L) _1//s vs. coupling strength J for 
the triangular lattice, with Xh = 1/8 . The lines connecting the data points are added for clarity. 

IV. PERCOLATION ANALYSIS 

For each spin configuration generated by the Monte Carlo algorithm, we performed a full 
decomposition in percolation clusters, using the randomly oriented Cartesian frame in the 
spin space as chosen in the preceding Monte Carlo step, and then placing bonds between 
nearest-neighbor pairs with probabilities Pij = max(0, 1 — e~ 2KSiS i ). The variable parameter 
K > governs the percolation process. While these percolation clusters are not involved 




S 



in spin-updating, they reduce to those obtained during the cluster simulations in the case 
K = J. To analyze this percolation problem, we sampled several quantities, including the 
second and fourth moments S 2 and S4 of the cluster size distribution, the Binder ratio Qi, 
the correlations g p (r = L/4), g p (r = L/2), and the ratio Q p . In this Section we describe the 
numerical results and analyses, and also an exact result for the present percolation problem 
on the triangular lattice. 



A. Percolation on the square lattice 

1. High-temperature range 

For tanhK = 1, i.e. in the limit K — > 00, all pairs of nearest-neighbor spins are connected 
as long as their x components are pointing in the same direction. At zero coupling J = 0, 
spins at different sites are uncorrelated, so that the percolation process reduces to standard 
site-percolation process, since the site occupation probability p — 1/2 may be identified 
with the sign of s x . An unimportant difference is that the present process forms percolation 
clusters for all the lattice sites while the standard site percolation constructs clusters only 
for the occupied sites. The site-percolation threshold p s c on the square lattice is very close 



to 0.592746 32n34|, and thus no infinite percolation cluster can occur at zero coupling 



strength J = 0, even for tanhA" = 1. Furthermore, from the results in Ref. [2lf| , where 
a similar percolation problem is studied in the context of several Potts models, we expect 
that no percolation transition occurs on the square lattice for small J. This expectation 
was confirmed by Monte Carlo simulations that were performed at several nonzero J < J c . 
Variation of K did not yield any signs of a percolation threshold. 



2. Low-temperature range 

The low-temperature XY phase J > J c displays algebraically decaying spin-spin corre- 
lations, which, unlike the exponential decay at J < J c , allows the formation of a divergent 
percolation cluster for sufficiently large K. We may thus expect a percolation threshold to 
occur at a J-dependent value K C (J). 

The existence of a percolation threshold K C (J) for J > J c is shown by the intersections 
of the curves in Fig. [6l which displays Q p as a function of K at J = 3.0 for several L. 
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FIG. 6: (color online). Dimensionless ratio Q p vs. the parameter K inducing the percolation 
transition. These data apply to the model on the square lattice, with spin coupling J = 3.0. The 
lower figure is an enlarged version. The lines connecting the data points are added for clarity. 



These data show that K C (J = 3.0) ~ 0.505. For K < K C (J), Q p rapidly approaches zero, 
as expected from the absence of long-range correlations of g p (r). 

In view of the long-range spin-spin correlations for J > J c , we have no reason to expect 
that the percolation transition at the threshold K C (J) for J > J c belongs to the uncorrelated 
percolation universality class. This is supported by the observation that, at J = J c , the 
fractal dimension of clusters with K = X is 2 — Xh = 15/8, which is different from the value 
91/48 for critical percolation clusters [351 ] . A closer look at the plot of Q p vs. K (Fig. [6]) 
indicates that, for K > K C (J), Q p rapidly converges to a independent nontrivial value 
smaller than 1. We propose the interpretation that, like the thermal transition induced by 
the variation of J, the percolation transition induced by K is also BKT-like. 

In Fig. [7] we display the correlation g p over a distance r = L/2 as a function of the linear 
system size L for several values of K. This figure shows a dependence of g p on L that 
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FIG. 7: (color online). Correlation g p {L/2) vs. linear system size L for several values of K which are 
shown in the inset. The use of logarithmic scales displays the approximate power-law dependence 
on L. These data apply to the square-lattice XY model at J = 3.0. The lines connecting the data 
points are added for clarity. 
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FIG. 8: (color online). Scaled correlation g , p L 1//4 (lnL) -1 / 8 over a distance r = L/2 vs. K for 
various linear system size L shown in the inset. These data apply to the square-lattice XY model 
at J = 3.0. The lines connecting the data points are added for clarity. 

approaches power-law behavior for large L. This suggests that percolation clusters remain 
critical for K > K c . Furthermore the exponent governing the scaling of g p (r) appears to 
depend on K. 

Next, we fitted the g p {r = L/2) data at J = 3.0 by 

g p = L- 2 ^(g + g l L- 1 + g 2 L- 2 ) , (9) 
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FIG. 9: (color online). Phase diagram of the square-lattice XY model in the J — K parameter 
space. The horizontal line represents the thermal BKT transition, and the diagonal line applies 
to K = J, where the percolation clusters are just those formed by the cluster algorithm. The line 
connecting the data points is added for clarity. 

where Xh is the associated scaling dimension. The terms with gi and #2 describe the finite- 
size corrections, and the correction exponents are simply set at —1 and —2, respectively. 
The results are shown in Table [B 

We conjecture that the fractal dimension 2 — Xh of the percolation clusters at K C (J) 
assumes the exact BKT value with xh = 1/8. This conjecture is based on the BKT-like 
behavior of the percolation transition in the low-temperature range, and on the numerical 
evidence for Xh obtained from the correlation g p (r = L/2) in Table [H First, the percolation 
in the low-temperature range seems to be BKT-like. Second, the fit results for the scaling 
dimension Xh, when interpolated to K c as given in Table Hit yield a value close to 1/8. The 
data for the scaled quantity g p (r = L/2)L 1 / 4 (lnL) _1/ ' 8 , shown in Fig. [H]for J = 3.0, confirm 
the existence of intersections, apparently converging to the same value of K as those in Fig. El 
Furthermore we found that the data for g p (r = L/2) in the interval 0.47 < K < 0.515 at 
J = 3.0 are well described by Eq. dSJ) for finite sizes in the range 32 < L < 400. This fit 
yields an estimate for the percolation threshold at K C (J = 3.0) = 0.504 (8). 

We also performed simulations at J = 2.4, 1.8, 1.4, and 1.2, and observe a behavior 
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FIG. 10: (color online). Scaling dimension Xh for percolation clusters at various points (J,K) for 
the square lattice. The values of J are shown in the inset. The horizontal line is located at Xh = 1/8 
and corresponds with the critical line K C {J). This figure indicates that the scaling dimension Xh 
depends on K as well as on J, and approaches 1/8 for K — > K C (J). The lines connecting the data 
points are added for clarity. 

similar as that described above for J = 3.0. On the basis of a fit of the (7 p L 1//4 (lnL) -1 / 8 data 
by Eq. (jSJ), we obtain the associated percolation thresholds, which are shown in Table HT1 
Next, we fitted Eq. (J9j) to the g p {r = L/2) data at various points (J,K) for J > J c and 
K > K C (J). The results are shown in Table HTT1 

In addition, we carried out simulations at J = 1.12, very close to the thermal critical 
point J c = 1.124 (3). The ratio Q P (K) appears to behave similarly as in Fig. [61 which 
suggests a BKT-like percolation transition. The estimated threshold K c = 1.120 (9) agrees 
with the critical point J c . This fits well with the continuation of the K C (J) line in Fig. [91 
Further, the numerical result for the fractal dimension of the percolation clusters at K = J c 
is consistent with the BKT value 2 — Xh = 15/8. 

The results in Table [Til and II III are summarized in Figs. [9] and [TO], respectively. 
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TABLE I: Results for the scaling dimension at J = 3.0 for various values of K for the square 
lattice. Parameters L m in and L max are the minimum and the maximum system size between which 
the Monte Carlo data of g p are included in the fit. 



K 


3.0 


1.8 


1.2 


0.9 


0.72 


-^min 


20 


20 


16 


16 


12 




400 


400 


400 


400 


400 


Xh 


0.02916 


(2) 0.0345 (6) 0.0412 (6) 0.0492 (6) 0.0590 (6) 


K 


0.60 


0.555 


0.525 


0.510 


0.495 


-^min 


16 


16 


16 


16 


16 




400 


400 


400 


400 


400 


Xh 


0.0751 (6) 0.0857 (6) 0.0990 (6) 0.1108 (6) 0.1290 (8) 



TABLE II: Percolation threshold K C {J) for various values of the XY coupling strength J. 



square 


J 

K c 


3.0 2.4 1.8 1.4 1.2 1.12 
0.504 (8) 0.553 (5) 0.646 (9) 0.785 (9) 0.946 (8) 1.120 (9) 


triangular 


J 

K c 


2.2 1.8 1.4 1.0 0.8 0.6824 
0.300 (3) 0.310 (1) 0.356 (4) 0.424 (9) 0.520 (9) 0.675 (9) 



B. Percolation on the triangular lattice 

1. Matching property 

The matching property 0, [37I plays an important role in the determination of the site 
percolation thresholds of several two-dimensional lattices; here we briefly review this subject. 
For a given planar lattice V = (V, B), where V is the set of lattice sites and B is the edge 
set, one does the following: 1), select parts of the faces of V, and fill in all the "diagonals" 
in those faces. This yields lattice C = (V, B + A), where A represents the set of all added 
diagonal edges. 2), select the faces that are not picked up in step 1), and fill in all the 
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TABLE III: Results for the scaling dimension Xh at various points (J, K) for the square lattice. 



J = 3.0 


K 

Xh 


30 18 12 09 72 
0.02916 (2) 0.0345 (6) 0.0412 (6) 0.0492 (6) 0.0590 (6) 


K 

Xh 


0.60 0.555 0.525 0.510 0.495 
0.0751 (6) 0.0857 (6) 0.0990 (6) 0.1108 (6) 0.1290 (8) 


J = 2.4 


K 

Xh 


2.4 1.44 1.2 0.96 0.84 
0.03748 (5) 0.0460 (6) 0.0491 (6) 0.0560 (6) 0.0628 (6) 


K 

Xh 


0.72 0.648 0.60 0.576 0.552 
0.0719 (6) 0.0821 (6) 0.0925 (6) 0.1010 (6) 0.1151 (6) 


J = 1.8 


K 

Xh 


1 8 1 26 1 08 9 81 
0.05282 (5) 0.0614 (6) 0.0672 (6) 0.0759 (6) 0.0832 (6) 


K 

Xh 


0.72 0.684 0.666 0.648 
0.0954 (6) 0.1028(6) 0.1069 (6) 0.1140 (6) 


J = 1.4 


K 

Xh 


1 4 1 19 98 854 798 
0.07386 (5) 0.0805 (6) 0.0915 (6) 0.1033 (6) 0.1125 (6) 


K 

Xh 


0.784 0.77 0.756 
0.1157 (6) 0.1202 (6) 0.1250 (6) 


J = 1.2 


K 

Xh 


1.8 1.44 1.2 1.14 1.08 
0.0798 (6) 0.0873 (6) 0.09610 (6) 0.0988 (6) 0.1028 (6) 


K 

Xh 


1.02 0.96 
0.1084 (6) 0.1159 (6) 


J = 1.12 


K 

Xh 


3.0 2.5 2.0 1.5 1.21 
0.0784 (3) 0.0821 (3) 0.0876 (4) 0.0984 (5) 0.1108 (3) 


K 

Xh 


1.18 1.15 
0.1132 (3) 0.1162 (5) 
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diagonals in these faces. One has lattice then C* = (V, B + A*) with A* the set of diagonals 
drawn in step 2). One calls lattices £ and C* are matching to each other; note that C and 
C* may be non-planar. Since no "diagonal" can be filled in a triangle, the triangular lattice 
is self- matching. It can be shown that, for the site-percolation problem, the cluster numbers 
per site n(p) and k*{1 — p) on a pair of matching lattices C and C* satisfy 

K (p)-K*(l-p) = <f>(p), (10) 

where p is the site-occupation probability and <fi(p) is a finite polynomial (it is termed 
"matching polynomial"). Equation (flOj) indicates that, if the cluster-number density k on 
the lattice C exhibits a singularity at a site occupation probability p, the same singularity 
will also occur in k* on C* at 1— p. Together with the plausible assumption that there is only 
one transition, the matching argument yields that the percolation threshold is p c = 1/2 for 
all self-matching lattices like the triangular lattice; further, it requires that <f)(p = 1/2) = 0, 
which is indeed satisfied by the result [37]] (f>(p) = p(l —p)(l — 2p) for self-matching lattices. 
An important feature of the matching argument is that it is still valid in the presence of 
interactions, as long as these interactions are symmetrical under the interchange of occupied 
and unoccupied sites. 



2. Percolation at tanh K = 1 

As mentioned in Sec. IIV A 1\ the case tanlii^ = 1, J = in the present percolation 
process corresponds with the case p = 1/2 for the standard-site percolation. The standard 
site-percolation threshold for the triangular lattice is p c = 1/2, thus the percolation threshold 
of the present percolation problem at J = is t&nhK = 1. 

Since the matching argument is independent of the coupling J, and no spontaneous 
symmetry breaking occurs in the two-dimensional XY model, tanh K = 1 describes a critical 
line for finite J. Further, we expect that, in the high-temperature range J < J c , the 
percolation transition is in the universality class of standard uncorrelated percolation, since 
there is no long-range spin-spin correlation. 

Figure [TT1 shows the data for the ratio Qi, which is defined by Eq. ([5]) on the basis of the 
size distributions of the percolation clusters. For J > J c , Qi approaches a J-dependent value 
which is clearly smaller than 1. This implies the absence of an infinite cluster that occupies 
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FIG. 11: (color online). Ratio Qi vs. coupling strength J for the triangular lattice at tanhA' = 1. 

The lines connecting the data points are added for clarity. 
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FIG. 12: (color online). Ratio Q p vs. coupling strength J on the triangular lattice at tanhX = 1. 
The lines connecting the data points are added for clarity. 
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FIG. 13: (color online). Scaled correlation g p (L/2)L 2xh vs. coupling strength J for the triangular 

lattice at tanh K = 1, with = 5/48 which applies to the uncorrelated percolation universality 

class. The lines connecting the data points are added for clarity. 

17 



a finite fraction of the whole lattice. The singularity at the thermal transition point J c is 
reflected by the jump that develops near J c w 0.68. 

Figure [TJ] shows the data for the ratio Q p . For J < J c , Q p converges to a universal 



value Q pc ~ 0.95 (note that this value differs from Q pc = 0.872776 (3) [38j for standard 
percolation on the triangular lattice, due to the difference mentioned in the first paragraph 
of Sec. IIV A II) . For J > J c , Q p approaches a J-dependent value smaller than 1; we thus 
expect that the correlation g p (r) decays algebraically rather than exponentially. The thermal 
transition at J c is reflected by the rapid variation of Q p near J c rj 0.68 for large L. 

The data for the scaled correlation g p (L/2) at tannic = 1 are shown in Fig. [T3] as a 
function of J, with Xh = 5/48 for the uncorrelated percolation universality. The convergent 
behavior for J < J c as a function of L confirms that the transition in this range belongs to 
the standard percolation universality class. The intersections roughly represent the thermal 
transition point J c . 



3. Percolation at tanh K ^ 1 

Following similar procedures as in Sec. IIV we obtain a percolation line K C (J) in low the 
temperature range J > J c for the triangular lattice. For J < J c we do not find a percolation 
threshold at finite values of K. The numerical results are shown in Table HT1 and Fig. [TH It is 
observed that the percolation in the range J > J c is also BKT-like, with a fractal dimension 
2 — Xh = 15/8 at K C (J), and a scaling dimension Xh depending on parameters K and J in 
the range K > K C (J). This is consistent with the results for the square lattice in Sec. IIV Al 



V. DISCUSSION 



Since spins in the same cluster formed during the simulations must have x-components of 
the same sign, the absence of a spontaneous magnetization [l| in the XY model means that 
the density of the largest cluster in the thermodynamic limit is also restricted to be zero, at 
least for finite values of J. The same restriction thus applies to percolation clusters formed 
with K = J in Eq. ([2]), and it must also hold for K < J. The absence of a nonzero density 
of the largest percolation cluster is in agreement with the interpretation of the percolation 
transitions for K < J described in Sec. IIVI as BKT-like. 
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FIG. 14: (color online). Phase diagram in the J — K parameter space for the XY model on the 
triangular lattice. The horizontal line is the thermal BKT transition, and the diagonal line is for 
K = J, where the percolation clusters are just those formed by the Monte Carlo cluster algorithm. 
In addition there is a percolation line for J < J c at tanhii" = 1, outside the range of this figure. 
The line connecting the data points is added for clarity. 



The results presented in Sec. IIV Al for the square lattice indicate that, in addition to the 
line K C (J) with J > K, also the line J = J c , K > J c represents a percolation threshold. As 
can be seen from the data points for J = 1.12 « J c in Fig.dUl the magnetic exponent depends 
on K along the latter line, which is thus a "nonuniversal" line of percolation transitions, 
and the bond dilution field parametrized by K is truly marginal. The existence of a BKT 
transition induced by varying K at the point K = J = J c corresponds with a marginally 
relevant bond-dilution field in the K < J c direction. 

As mentioned earlier, for the triangular lattice with J < J c , the line tanhi^ = 1 is 
critical, and belongs to the standard percolation universality class. The continuation of the 
tanh K — 1 line to J > J c is also critical (in the sense that the correlation functions display 
algebraic decay), but with a J-dependent critical exponent. 

In order to obtain some more information on the dependence of the present percolation 
problem on the coordination number z, we also simulated the z — 18 equivalent-neighbor 
XY model on the triangular lattice, which has equal nearest-, second-nearest- and third- 
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nearest-neighbor interactions. The procedure outlined in Sec. (JTTTJ) yielded an estimate of 
the thermal transition at J c = 0.162 (2). As expected, a Monte Carlo analysis of the 
percolation problem with z = 18 showed the existence of a critical line K C (J) in the high- 
temperature phase J < J c , belonging to the standard percolation universality class. When 
J approaches J c , the K C (J) line bends toward large values of K. This suggests that K C (J) 
line ends at K — > oo for J = J c . For K > K C (J) there is clear evidence for the existence of 
a percolation cluster with a finite density in the limit of large L. 

In the low-temperature range J > J c , we found, just as for the models with nearest- 
neighbor interactions, a BKT-like transition line K C (J), as in Figs. IHlandHH In spite of the 
relatively large coordination number z = 18, no evidence is found for a percolation cluster of 
a nonzero density, even at tanhfT = 1. Although the spin-spin correlations in the algebraic 
XY phase for J > J c stimulate the percolation transition in the sense that it occurs at 
smaller values of K when J increases, it also appears that they obstruct the formation of a 
percolation cluster with a non-zero density. 

Besides the rule based on Eq. (j2J), other procedures for placing bonds may be applied. 
For instance, as mentioned in the Introduction, Wang et al. placed percolation bonds 
between neighboring XY spins if their orientations differ less than a given threshold, and 
found percolation transitions in the uncorrelated percolation universality class for all XY 
couplings. Another possibility is to place percolation bonds with probabilities given, instead 
of Eq. ([2]), by pij = max 



similar to those of Ref. 
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(0,1 - e - 2Kgi "'). In that case, we expect percolation transitions 



including transitions in the low-temperature range J > J c of the 



XY model. Indeed, a preliminary Monte Carlo analysis of this problem 2 



confirms the 



existence of such transitions in the universality class of uncorrelated percolation. 

Finally we remark that recently a percolation problem was formulated on the basis of 
the O(n) loop configurations 22j]. Since the XY model is equivalent with the 0(2) model, 
another way thus arises to introduce percolation in XY-type models. Although this seems 
a very different approach, it reproduces our result that a marginally relevant dilution field 
exists at the BKT transition. 
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